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Abstract 

In this paper we give definitions of matrix rates of return which do not depend on the choice 
of basis describing baskets. We give their economic interpretation. The matrix rate of return 
describes baskets of arbitrary type and extends portfolio analysis to the complex variable domain. 
This allows us for simultaneous analysis of evolution of baskets parameterized by complex variables 
in both continuous and discrete time models. 
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I. INTRODUCTION 



The goal of capital investment is the maximization of a profit and minimization of 
possible losses. This goal cannot be achieved by investing of the whole capital in the most 
profitable enterprises. Such situations do not happen. The future profit of a market in- 
vestment is uncertain, therefore, the investor creates composite baskets consisting of capital 
investments of possibly most diversified character. This kind of procedure diversifies the 
risk of enterprise. The description of the evolution of multidimensional capital of this kind 
is essential to quantitative analysis of the correlations related with investment processes, in 
particular, these for which the tools of traditional financial mathematics are unapplicable. 
Depending on the point of view, in every capital basket we can, besides of quantitative 
changes of individual components, observe the flows between its components. The flows of 
capital can be recorded even if no decisions about capital operations are taken. Such situ- 
ations require matrix description and enforce the generalizations of the calculus of interest 
rates in such a way that they are sensitive to both quantitative changes of individual ele- 



ments and flows between them, 
of the interest rate calculus, see 



his situation leads inevitably to the matrix generalization 
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II. LINEAR HOMOGENEOUS CAPITAL PROCESS 

Let us consider a capital of a banker which is a mixture of two elements: k\ - an amount 
lent to client, k 2 - the remaining of assets. Let the variable I G Z denote arbitrary time 
interval. The process of change of the banker's capital related with the amounts k\ and k 2 
can be described as follows 

ki(l): According to the rate of interest «i(7) at which the banker gave credit, the growth of 
the component k\ is equal to ai(l)k\(l). In addition the amount of unreturned part of 
credit decreases by the value of repayment (3{l)ki{l) which is determined by the rate 

k 2 (l): The capital k 2 increases by the amount of repayment of credit f3(l)ki(l). Being in 
addition placed in, for example liquid stocks of annual rate of return a 2 (l), it grows 
by the amount a 2 (l)k 2 (l). 
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The components of the banker's capital do form a basket k(l) = (k\(l), k 2 (l))- It is repre- 
sented by an element of the two dimensional real vector space R 2 . The evolution of this 
basket during the time I G Z can be described by the system of equations 

k\{l + 1) = (1 + «i(Z) - (3(1)) ki(l) 
k 2 (l + l) = (3(1) k x (l) + (1 + a 2 (l)) h(l) ' 

We interpret the negative values of components of the basket k m , m = 1, 2, . . . as debts of 
the banker. 

The nonlinear (with respect to the remaining debt) repayment rates, that is for example 
the costs of service of credit or taxes, can be presented in a form of linear repayment, after 
appropriate modification of the factor f3(l). In particular various borrower's obligations to 
the banker can be taken into account with the help of the expression in variable proportion to 
the amount of paid-off credit. This enables us to apply our formalism in much wider context. 
This evolution of the capital can be described in such a way that any changes of components 
of the basket are expressed as percent changes of those components. In the vector space of 
baskets we can choose a new basis such that we do not observe any flows of capital, but 
autonomous growth of individual components only. Let us assume that a.\(l) ^ a 2 (l) and 
k x (l) is as in previous basis, the amount of the capital lent to client k[(l) = k\(l) . k' 2 (l) is 
the sum of (3(1) part of a loan and ot 2 (l) — ol\(1) + (3(1) part of remaining banker's assets. 
That is k' 2 (l) = (3(l)ki(l) + (a 2 (l) — ol\{1) + (3(l))k 2 (l). In these new variables equations (0), 
describing the changes in the basket separate and take the following form: 

k[(l + l) = (l + ai(Z) -(3(1)) k[(l) 
k' 2 (l + l)= (l + a 2 (l))k' 2 (l) ' 

In this way we do not observe any flows of capital and the components k\ and k' 2 of the 
basket grow according to the rates of interest a±(l) — (3(1) and a 2 (I) for k[ and k' 2 respectively. 
The formalism of financial mathematics should not depend on choice of basis to describe the 
baskets. It seems that the matrix rates presented below are necessary to obtain the basis 
independent description of the economical reality. 
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III. THE MATRIX RATE OF RETURN 



Capital processes are described by linear homogeneous difference equations. In the matrix 
notation they take the form: 

k(Z + l) -k(Z) = R(Z)k(Z), hence k(Z + 1) = (I + R(Z)) k(Z) , (2) 

where k(Z) e K M , R(Z) and I are real matrices of dimension M x M. I denotes the unit 
matrix. Matrix R(Z) is called the matrix credit rate of return or the matrix lower rate. The 
adjective lower refers to the way of the description of the growth of the vector k(Z). These 
changes are expressed as the effects of linear transformations applied to the vector k(Z) of 
the basket in the moment Z prior to this change. The space R M of all possible baskets, is 
called the phase space of baskets and linear homogeneous difference or differential systems 
of equations of first order are called equations of motion of the basket. For M = 2 the matrix 
R(Z) generating the evolution ((TJ) equals to 

\ W) Mt) ) 

If the state of the basket in the initial moment p equals to k(p), then the solution of the 
evolution equation (J2J) takes the form 

k(r)= (rn(I + E(«)))k(p), (4) 

where T denotes the chronological ordering operator. It orders all matrices chronologi- 
cally: the matrices with later arguments are moved to the left of the matrices with earlier 
arguments, that is 

p+n 

T Y[ A(s) = A(p + n)A(p + (n - 1)) . . . A(p + l)A(p) , 

s=p 

where A(s) is an arbitrary sequence of matrices. 

The chronological product T[]^~ (I + R(s)) contains more detailed information about 
changes of capital in the basket than the usually used quotients ^JT 1 fc m ^ . The possibility 
to use the standard calculus of interest rates depends on way of description of the basket, 
that is on choice of a basis in our language. For this reason, in order to obtain an invari- 
ant description, the rate of growth R(Z) is promoted to a composite object with definite 
transformation rules corresponding to the changes of reference frame (an observer). 



Let us consider processes for which the matrix I + R(Z) is nonsingular. Then, introducing 
the concept of the matrix discount rate of return or the matrix upper rate R(Z), it is possible 
rewrite the equation of motion (j2J in the following form: 

k(Z + 1) - k(Z) = R(Z) k(Z + 1) . (5) 

Comparing formulae (J2J) and (0) we obtain the relation between both types of matrix rates 
introduced above: 

(I + R(Q) (I - R(Z)) = (I - R(Z)) (I + R(Z)) = I (6) 

Solving the above equation with respect to R(Z) we obtain 

R(Z) = R(Z)+R 2 (Z)+R 3 (Z) + .... (7) 

From the point of view of capitalisation from upper, we can interpret the above formulae as 
the contribution of discount rate into increase of capital by summation of all interests from 
interests (the geometrical sequence). According to Eq. (jZJ) the formulae for R(Z) takes the 
form -R(Z) = (-R(Z)) + (-R(Z)) 2 + (-R(Z)) 3 + . . . . Notice that for fixed argument I the 
matrix rates are commuting. 

When we perform the formal change of the direction of time the credit and discount rate 
change their signs and they change their roles. Therefore, the formulae containing these 
matrices are symmetric with respect to the time reflection. The appropriate formulae for 
(HJ) can be obtained with the help of (JHJ) as the solution of the equation of motion: 

k(p)= (r'n(i-R>))W) ; (8) 

where T' is the antichronological operator which orders matrix rates in the direction opposite 
to that corresponding to the operator T . 

If the matrix rates at different moments of time commute, that is in the case of M = 1 or 
for any M if matrix rates are time- independent, we can neglect the operators T and T in 
the solutions and (jHJ. 

IV. INTERPRETATION THE MATRIX RATE OF RETURN 

The matrix rate of return R(Z) = (R mn (Z)) can be given as sum of two matrices R(Z) = 
C(Z) + D(Z), with C(Z) the matrix of flows (the name is justified by the property that sum 
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of the elements of each column is equal zero) and D(Z) - the diagonal matrix called the 
matrix of growths. This decomposition is unique, when the basis is fixed. Introduction of 
the matrix rate of return is essential, when we cannot transform the matrix C(/) to matrix 
zero. Moreover, the off-diagonal elements of the matrix R(Z), R mn (0 for m ^ n, see Eq. 
(jHJ, describe which part of capital of the n-th component of the basket flows to the m-th 
component of the basket. Diagonal elements R mm (Z) describe the growths. The rate of 
growth of the m-th component of the basket equals a m (l) = R mm (Z) — J2 m ^n^nm^)y that 
is, the diagonal element corrected by all out-flows of the capitals related to the component 

MO- 



A. Example 

(i) Let us restrict the process from section II to the cases when 0(1) = and the coefficients 
ctj do not depend on time I and the second coefficient is two times bigger than the first 

, a 

one. Then the matrix lower rate is equal to R = . According to the principal 



2a 



formula © the corresponding matrix upper rate is given by: R 



1+Q 







2a 



l+2a 

The process contains the matrix of growths only because the matrix of flows is zero. 
It is possible to analyze the profits on the grounds of autonomic evolution in one- 
dimensional subspaces of the phase space - the classical concept of interest rates is 
applicable here. 

(ii) We can look at this process in a different way. Namely by describing it in the coordinates 
in the new basis in the space of baskets. Let the reference basis consist of client's debts 
and all banker's capital assets. Then the equation of motion of the basket (JIJ takes 
the form 

'k 1 (i + l) = (l + a)k 1 (Z) 
k 2 (/ + l)= -aki(0 + (l + 2a)k 2 (/). 
The matrix of flows determined by the above equation is non zero now. The matrix 
of flows and the matrix of growths equals to 

C 



/ a 1 




^0 







I- H 






1 -a 1 




\o 


2a j 



respectively. The debt of a client k\ changes in a similar way to the previous example, 
though not as result of an autonomic growth, but due to an outflow of the banker's 
capital. Banker's capital grows according to the same rate as investment k 2 (l)- 

(iii) In the basis given by Eq. (JHJ) the matrix upper rate for the equation of motion (JHJ) 
takes the form 



Different points of view at the same capital process presented in (i), (ii), (ii) are equally 
correct and sensible. Comparing the credit convention in variant (ii) with the discount 
convention (iii) we can note the essential difference between the matrices D and D. Growth 
of first component is only the effect of flows in the first case while in second case this 
component has partial autonomy in its growth. The indicated difference in interpretation 
can be a reason of many financial embezzlement, exactly in the same way as it happens in 
case of nonpayment of interests of overdue interests in simple capitalisation. We will refer 
the asymmetry of this kind in the description of flows as the paradox of difference rates. 

The interpretation of the matrix rate R is analogous to that of the matrix R with the 
only difference that the rate R defines the capital changes with respect to the oncoming 
moment of the settlement of accounts. The convention of discount replaces the convention 
of credit. The free choice of convention is not possible for nonreversible processes i.e. when 
the matrices R and R are singular. 

V. THE FORMALISM OF CONTINUOUS DESCRIPTION OF CREDIT 

For the capital calculus to be transparent and readable for practitioners the formal so- 
lution of the equation of motion (j2J) can be modeled numerically or presented it in dense 
form with the help of the limiting process which transforms the discrete models described 
by linear difference equations into continuous ones with differential equations. Assume that 
we consider time scales such that the periods of time between the changes of components of 




and it is the sum of the following matrices of flows and growths 
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the basket is infinitesimal and equals to r = ti+i — U. After rescaling of the time domain of 
the basket the equation of motion (J2J) takes the form: 

k(t, + r) - k(t t ) R(Z) 



k(t,) . 



In the limit r->0we obtain 



= R(t) k(t) where R(t) := lim ^-\ t =t r (10) 

dt r^O T 

Matrix R(£) is called the differential matrix rate of return. 

The formal solution of the equation (fTUj) is given by the formula 

k{t) = (W'o R(/) *'Wi ). (ii) 

The chronologically ordered exponential function is infinite series in the differential matrix 
rate 

Te f to m)dt =1+ / R(t 1 )dt 1 + / R(^) / R(t 2 )dt 2 dt 1 + .... 

J to J to J to 

If the differential matrix rate is constant, the chronological operator T is the identity The 
expression on the right hand side of the equation (fTTj) which describes the time evolution of 
the basket becomes transformed to the standard matrix exponential function Te^o R( '* ^ dt = 

g(t— to)Ro 



VI. COMPLEX RATE OF RETURN 



Let us consider baskets which have oscillating components. Part of their capital becomes 
a debt at some moment and a desirable good at another time. This phenomenon is called 
pumping of capital in the language of financial market. Evolution of basket of this kind is 
described by one of the most popular model of physics - the harmonic oscillator. Let us 
consider the differential matrix rate 

\b o J yo -bj \b b J 

It is easy to imagine the corresponding bank contracts leading to the flows of capital in basket 
and autonomic increases of components of the basket generated by these flows. Let us con- 
sider the complex extension C 2 of the phase space M. 2 . Then kx = (1, i) and k 2 = (1, — i) = k* 
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are the eigenvectors of the matrix rate R(t) with eigenvalues — i b and i b respectively. The 
description of process is simplified because the basket splits into two independent compo- 
nents with abstract complex capital. Although that baskets of the complex capital are the 
abstract concepts, they have the interpretation in the real basis due to matrix of transition. 
The equation of motion of the basket in the basis of eigenvectors {k 1; k 2 } takes the form: 

k>) = e- ib ^ kx(to) , k 2 (t) = e^-^ k 2 (t ) . 

In the initial, real basis one has 

/ cos(6(*-*o)) -sm(b(t-t ))\ 
k(t) = k(t ) • 

\^ sin(6(t - to)) cos(6(*-*o)) / 

This equation describes motion along the circle centered at the beginning of the cartesian 
coordinates of the basket. The period of return to the same point of the phase space equals 

to T = m- 



VII. THE INDEFINITE MATRIX RATE 

We define the indefinite logarithm, see 0], to be a mathematical object representing the 
abstract concept of the logarithm with an unfixed base. For any given real number x>0, the 
indefinite logarithm of x written as [log re], is a special type of mathematical object called a 
logarithmic quantity object, which we define as follows [2]: 

[log a;] := {(b,y)\b > 0,y = \og b x} . 

Indefinite logarithmic quantities are inherently scale-free objects, that is, they are non-scalar 
quantities and they can serve as a basis for logarithmic spaces, which are natural systems 
of logarithmic units suitable for measuring any quantity defined on a logarithmic scale. 
Although the above definition is restricted to positive real numbers, it could be extended 
to non-zero complex numbers too. The concept of the rate of interest is connected with the 
time scale. To get rid of explicit time scale dependence we introduce the indefinite matrix 
rate, which is the generalization of the indefinite logarithm to the multidimensional case: 

togfVe^V'y . 

In the above definition the argument of the logarithm is a matrix. 
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VIII. CONCLUSIONS 



Every matrix can be deformed to a diagonalizable complex matrix by arbitrarily small 



deformations of their elements 



4 1 , that is the set of diagonalizable linear transformations 
of the complexified phase space C M is dense in the space all linear maps of C M . From the 
above it follows that, the evolution of every capital basket can be represented as set of non 
interacting complex capital investments. Therefore, in the complex extended phase space the 
decomposition of the matrix rate R(t) = C(t) +D(t) can always be done in such a way that 
the matrix of flows is zero. As in traditional financial mathematics the real parts of nonzero 
elements of a diagonal matrix rate measure the loss or gain of complex investments. The 
imaginary parts inform about periodicity of changes in proportion between real component 
of complex investments and its imaginary partner. To every periodic changes of proportion 
of the elements of the complex components of capital there will correspond its complex 
conjugated partner. The evolution of real capital of basket is most easily observed in terms 
of its components with respect to the basis of eigenvectors of matrix rates. 

By the choice of appropriate moments of entering and exit from capital process the 
oscillations like above can be used as particularly effective mechanism of enlarging of the 
capital giving the similar results as financial leverage. 
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